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Many papers have been written on theories described by non-Hermitian VT- 
symmetric quantum-mechanical Hamiltonians. To construct such theories one begins 
with a Hamiltonian that is both Hermitian and VT symmetric, such as the harmonic 
oscillator H = p 2 +x 2 . One then introduces a real parameter e to extend the Hamiltonian 
into the complex domain in such a way as to preserve the VT symmetry: 

H = p 2 + x 2 (ix) e . (1) 

The result is a family of complex non-Hermitian Hamiltonians that for positive e 
maintain many of the properties of the harmonic oscillator Hamiltonian; namely, that 
the eigenvalues remain real, positive, and discrete [UI21IS]- The properties of classical 
PT-symmetric Hamiltonians have also been examined jUEllEllZI- However, there are 
to date no published studies of PT-symmetric classical wave equations. 

The starting point of this paper is the heretofore unnoticed property that the 
Korteweg-de Vries (KdV) equation, 

u t + uu x + u xxx = 0, (2) 
1 Permanent address: Department of Physics, Washington University, St. Louis MO 63130, USA 
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is VT symmetric. To demonstrate this, we define parity reflection V by x — > —x, and 
since u = u(x,t) is a velocity, the sign of u also changes under V: u —>■ —u. We define 
time reversal T by t — > —t, and again, since u is a velocity, the sign of u also changes 
under T: w — > — it. Following the quantum-mechanical formalism, we also require that 
i —>■ —i under time reversal. It is clear that the KdV equation is not symmetric under V 
or T separately, but it is symmetric under combined VT reflection. The KdV equation 
is a special case of the Camassa-Holm equation [S], which is also VT symmetric. Other 
nonlinear wave equations such as the generalized KdV equation u t + u k u x + u xxx = 
and the Sine-Gordon equation u tt — u xx + g sin u = are VT symmetric as well. 

The striking observation that there are many nonlinear wave equations possessing 
VT symmetry suggests that one can generate many families of new complex nonlinear 
PT-symmetric wave equations by following the same procedure that was used in 
quantum mechanics [see (|T|) ]. One should then try to discover which properties of 
the original wave equations are preserved. 2 

In this brief note we limit our discussion to the complex PT-symmetric extension 
of the KdV equation: 

Ut - iu(iu x ) e + u xxx = 0, (3) 

where e is a real parameter. We now examine the remarkable properties of some members 
of this family of complex PT-symmetric equations. We will emphasize the properties 
that members of this class have in common. 



Case e = 1: When e = 1, Q reduces to the KdV equation (J2J. The KdV equation has 
an infinite number of conserved quantities [9 . The first two are the momentum P, 

^P = 0, P = J dxu(x,t), (4) 
and the energy E, 

-^£ = 0, E = \ J ' dx[u(x,t)} 2 . (5) 

The Cauchy initial-value problem for the KdV equation can be solved exactly because 
the system is integrable, and it is solved by using the method of inverse scattering [TO] . 
A solitary-wave solution to the KdV equation has the form 

u(x, t) = 3csech 2 [|\/c(x — ct — x )] , (6) 

where c > is the velocity. (In general, a solitary-wave solution u(x, t) = f(x — ct) to a 
partial differential equation is defined to be a wave that propagates at constant velocity 
c and whose shape does not change in time. In this paper we also require that f(z) —>■ 
as \z\ — > oo.) These solitary waves are called solitons because as they evolve according 
to the KdV equation they retain their shape when they undergo collisions with other 
solitary waves. One can observe numerically how a soliton emerges from a pulse-like 

2 An alternative possibility for study is to examine inverse scattering problems and isospectral flow 
using PT-symmetric potentials u(x, t). We reserve this research direction for a future paper. 
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initial condition. For example, for the initial condition u{x, 0) = 3sech(x), we see in 
Fig. ^ that the pulse sheds a stream of wave-like radiation that travels to the left and 
gives birth to a right-moving soliton of the form in (JHJ). 




■0.5 



Figure 1. Birth of a soliton. As the initial condition u(x,Q) = 3sech(x) evolves 
in time according to the KdV equation, a soliton of the form @) moves to the right 
leaving a trail of residual radiation that propagates to the left. The function u(x,T) 
is plotted for T = 0, 0.8, 3.5, 7, and 14. 



Case e = 0: Setting e = in (JHJ) gives the linear equation 

u t -iu + u xxx = 0. (7) 

To solve the initial- value problem for this equation, we substitute u(x, t) = e lt v(x, t) and 
reduce it to Vt + v xxx = 0. We then perform a Fourier transform to obtain the solution 
in the form of a convolution of the initial condition and an inverse Fourier transform: 

v(x,t) =v(x,Q)*F- 1 (e ipH ^j . (8) 

The inverse Fourier transform of the exponential of a cubic is an Airy function [11J. 
Thus, the exact solution for u(x, t) is 

/oo 
dsu(x-s,0)Ai [(3t)~ 1/3 s] . (9) 
-oo 

The Airy function Ai(x) has a global maximum near x = 0. For x large and positive 
Ai(x) decays exponentially, Ai(x) ~ (2 v /7r)~ 1 x~ 1 / 4 exp (— |x 3 / 2 ), and for x large and 
negative Ai(x) decays algebraically and oscillates, Ai(— x) ~ n^^x" 1 ^ 4 sin (|^ 3 ^ 2 + ■ 
Thus, the qualitative behavior of Ai(x) resembles that in Fig. [T] If we choose the initial 
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condition u{x, 0) = 3sech(x) that was used to generate Fig. Q then apart from the 
phase e 1 *, we find that the solution ©, which is shown in Fig. El resembles that for 
the KdV equation in Fig. ^ except that no soliton emerges from the initial condition. 
There is only residual radiation that travels to the left. Thus we have extended the KdV 
equation into the complex domain while preserving many of its qualitative features. 
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Figure 2. Plot of u(x, T)e^ T as a function of x for T = 0, 5, 10, 20, 40, and 80, where 
u(x,t) in © is calculated for the initial condition u(x,0) = 3sech(x). A comparison 
of this figure with Fig. ^ shows that the solutions to 10 with e = and e = 1 are 
somewhat similar. However, in this case, while the initial condition produces radiation 
that travels to the left, it does not give rise to a soliton. 



Case e = 3: If we set e = 3 in Q, we obtain the nonlinear wave equation 

u t - u(u x f + u xxx = 0, (10) 

which has received only passing mention in the literature ^21 Ej • The only observations 
that have been made regarding this equation are that, apart from translation invariance 
in x and t, there is an obvious scaling solution of the form u(x, t) = <p (xt~ 1//3 ). Yet, this 
equation has an array of rich and beautiful properties that have so far been overlooked. 
As we will show, there are two conserved quantities, a momentum P and an energy E, 
which are the analogs of P and E in and (JSJ for the KdV equation. We will also 
show that there are traveling waves, and we will see how a pulse-like initial condition 
gives birth to a traveling wave, just as in the case of the KdV equation. 

We derive the conserved quantities for the wave equation (|1U|) in much the same 
way that one finds the conserved quantities for the KdV equation (j2J. However, the 
procedure is more elaborate. To find the momentum P, we begin by integrating (|rUjl 
with respect to x and assume that u(x,t) vanishes rapidly as \x\ — > oo. For the case of 
the KdV equation, this procedure immediately gives the result in (J3J). However, for the 
equation in (|10|) we have the result 

/ dxu = [ dxu(u x ) 3 . (11) 
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Evidently, f u is not a conserved quantity because the right side of of this equation 
does not vanish (in contrast to the KdV equation). To proceed we introduce the identity 

dxu N {u x f = — — — — / dxu N+2 u xxx , (12) 



(N+l)(N + 2) 

which is obtained by performing two integrations by parts. Using this identity for the 
case N = 1, we rewrite (fTTj) as 

d f 1 f 

dxu = — / dxu 3 u xxx . (13) 



dt J 3 _ 

This equation suggests that we should multiply (fTUjl by u 3 , integrate with respect to x, 
and use the identity (pH2j) for the case N = 4 to obtain 

- — f dxu A = — f dxu 6 u xxx — f dxu 3 u xxx . (14) 



Adt J 15 . 

We can combine (fT3*|) and fPfj) to eliminate the J u 3 u xxx term, but the right side will 
still not vanish because there will be a J u 6 u xxx term. 

Therefore, we must iterate this process by multiplying by u 6 , u 9 , u 12 , and so on, and 
then integrating with respect to x. We thus obtain the following sequence of equations: 

Jt\ dX WVl = \ dX (3kl 2)(3k X + 3) " / dX A — (15) 

where k = 0, 1, 2, 3, • • •. We can now completely eliminate the right side if we multiply 
the kth equation in (|T3j) by 

where A is an arbitrary constant, and sum from k = to oo. We conclude that = 0, 
where the conserved quantity P is given by 

By a similar argument, we can construct a second conserved quantity E, ^E = 0, 
where E is given by 

E = B / dxV ) - 3 \, 18 

J ^ (3^ + 2)! 

and B is an arbitrary constant. 

The summations in ()17j) and (fTSj) can be performed in closed form in terms of Airy 
functions, giving 

f2 x / 3 u{x,t) 

P = I dx I ds 

-2 1 / 3 n(a;,t) 



E 



Bi(s) + v / 3Ai(s) 

/p-^'"u(x,t) 
dx J ds [Bi(s) - VsAi(s) 



(19) 
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where we have taken A = 6 1 ^ 3 /n and B = 6 2 ^ 3 /n. It is especially noteworthy that 
the conserved quantity E is strictly positive when u(x,t) is not identically 0, and thus 
it is reasonable to interpret E as an energy. The positivity property of the energy is 
maintained when e changes from 1 (the KdV equation) to 3. We do not believe that 
(fTUj) has more than two conserved quantities. 

Equation (fTUj) is also similar to the KdV equation in that it has solitary-wave 
solutions. To construct such a solution, we substitute u(x,t) = f(x — ct) into (fTUj) to 
find the ordinary differential equation satisfied by f(z): 

-cf'(z)-f(z)[f(z)? + f"(z) = 0. (20) 

It is only possible to solve this autonomous equation in implicit form. To do so, we seek 
a solution of the form f'(z) = G(f). The function G satisfies 

-2c-2/G 2 (/)+[G 2 (/)]" = 0. (21) 

Making the further substitution H(f) = G 2 (f), we find that H satisfies 

H"(f) - 2fH(f) = 2c, (22) 

which is the inhomogeneous Airy equation, whose solution is expressed in terms of the 
inhomogeneous Airy or Scorer function Hi 11 . 

Unfortunately, because the solution to (j2*Uj) is implicit, it is not easy to determine 
immediately whether there are solitary-wave solutions [solutions f(z) that vanish as 
\z\ — ■> oo]. However, numerical analysis confirms that there are indeed such solutions. 
In Fig. El we have plotted the solitary wave for c = 1. Note that this wave is an even 
function of z and it decays like e~' z for large \z\. Note also that the solitary-wave 
solution is a negative pulse, rather than a positive pulse as with the KdV equation. 




Figure 3. Solitary-wave solution to the differential equation u t — u(u x ) 3 + u xxx = 0. 
This negative-pulse solution has the form u{x,t) = f(x — ct), where we have taken 
c = 1. The solitary wave is an even function of z — x — t and it decays like e _ ' z ' 
for large \z\. Thus, it closely resembles the solitary- wave solution in JJJJ for the KdV 
equation. At the negative peak the height of this solitary wave is —2.73802. 

As we saw in Fig. Q for the KdV equation, an initial pulse such as u(x, 0) = 
— 3sech(x) for ()10|) gives birth to a solitary wave. As shown in Fig. HI this initial pulse 
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emits radiation that travels to the left and evolves into a right-going solitary wave (see 
Fig. |3J). Computer experiments suggest that these solitary waves are not solitons; that 
is, they do not maintain their shape after a collision with another solitary wave. Indeed, 
we would be surprised if (II Oj) were an integrable system. The quantum-mechanical 
Hamiltonian(£Q) ceases to be exactly solvable when e ^ 0, and in the same vein we 
expect that © is not integrable when e ^ 0. 
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Figure 4. Birth of a solitary wave. As the initial condition u(a;,0) = — 3sech(x) 
evolves in time according to u t — u(u x ) 3 + = 0, a pulse that approaches the shape 
of a solitary wave moves to the right leaving a trail of residual radiation that propagates 
to the left. The wave u(x,T) is shown for the times T = 0.05, 0.25, 0.55, 1, and 2. 



There are no positive solitary- wave solutions to (jlOj) . As we see in Fig. an initial 
pulse of the form u(x, 0) = 3 sech(x) generates a stream of radiation that travels to the 
left, but it does not give rise to a solitary wave. 

Case e — 2n + 1: When e = In + 1 is an odd integer, the nonlinear wave equation in (j3J) 
is real: 

ut + (-l)M0 2n+1 + = 0. (23) 
For all values of n there are solitary waves u(x,t) = f(z), where z = x — ct, and these 
waves are even functions of z. As n increases, the solitary waves alternate between being 
strictly positive and strictly negative functions and gradually become wider (see Fig.EJ). 

In conclusion, we have shown how to extend the conventional KdV equation into 
the complex domain while preserving VT symmetry. The result is a large and rich class 
of nonlinear wave equations that share many of the properties of the KdV equation. In 
particular, we find that for some values of e there are conservation laws and solitary 
waves, and that arbitrary initial pulses can evolve into solitary waves after they give 
off a stream of radiation. Airy functions appear repeatedly in the analysis because the 
associated linear equation that is satisfied when u(x,t) is small is of Airy type. 
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Figure 5. The solution to u t — u(u x ) 3 +u xxx = that evolves from the initial condition 
u(x,0) = 3sech(x). This initial condition generates radiation that travels to the left, 
but it does not give rise to a solitary wave. The wave u(x, T) is shown for the times 
T = 0, 6, 12, and 19. 
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Figure 6. Solitary- wave solutions of the form f(x — t) to the differential equation 
u t + (— l) n u(u x ) 2n+1 +u xxx = for n = 1, 2, 3, 4. Note that the solutions are alternately 
positive and negative and gradually get wider as n increases. At the stationary points 
the heights of the waves are -2.73802 (n = 1), 2.45839 (n = 2), -2.30305 (n = 3), 
and 2.20797 (n = 4). The widths of the waves at half-height are 3.15, 3.14, 3.19, and 
3.26, respectively. 



In this paper we have just scratched the surface; one can begin with other VT- 
symmetric nonlinear wave equations, such as the Camassa-Holm or the generalized 
KdV equations, and study the properties of the resulting new complex wave equations. 
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